%% Carbon Taxes Around the World: Cooperation, Strategic Interactions, and Spillovers
% IMF Economic Review
% Alessandro Moro and Valerio Nispi Landi
% Replication files
% This file saves the IRFs under the Nash equilibrium

dynare climate
load path0
T=200;

variables_pct = {'gdp', 'gdpz', 'c', 'cz', 's', 'kB', 'yBz', 'kG', 'kB_zu', 'kG_zu', 'kB_zz', 'kG_zz', 'kB_zd', 'kG_zd'};
variables_dev = {'bY', 'tbY', 'ErkB', 'ErkG', 'ErkB_zu', 'ErkG_zu', 'ErkB_zz', 'ErkG_zz', 'ErkB_zd', 'ErkG_zd', 'r', 'rz', 'uip', 'mu', 'muz'};
variables_div100 = {'eW', 'e', 'ez', 'x'};

for i = 1:length(variables_pct)
    eval([variables_pct{i}, '2 = (', variables_pct{i}, 's(1:T))./(', variables_pct{i}, '0(1:T))-1;']);
end

for i = 1:length(variables_dev)
    eval([variables_dev{i}, '2 = ', variables_dev{i}, 's(1:T) - ', variables_dev{i}, '0(1:T);']);
end

for i = 1:length(variables_div100)
    eval([variables_div100{i}, '2 = (', variables_div100{i}, 's(1:T)-', variables_div100{i}, '0(1:T))/100;']);
end


save path2 gdp2 gdpz2 c2 cz2 kB2 kG2 kB_zu2 kG_zu2 kB_zz2 kG_zz2 kB_zd2 kG_zd2 ErkB2...
ez2 e2 ErkG2 ErkB_zu2 ErkG_zu2 ErkB_zz2 ErkG_zz2 ErkB_zd2 ErkG_zd2 s2 bY2 tbY2 eW2 x2 r2 rz2 uip2 yBz2 mu2 muz2
 
 